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In this paper, the single-spin transition dynamics is used to investigate the kinetic Gaussian 
model in a periodic external field. We first derive the fundamental dynamic equations, and then treat 
an isotropic d- dimensional hypercubic lattice Gaussian spin system with Fourier's transformation 
method. We obtain exactly the local magnetization and the equal-time pair correlation function. 
£h ■ The critical characteristics of the dynamical relaxation r q , the complex susceptibility x{ u} i c I)j and 

the dynamical response are discussed. The results show that the time evolution of the dynamical 
quantities and the dynamical responses of the system strongly depend on the frequency and the 
0^ ■ wave vector of the external field. 
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I. INTRODUCTION 

The purpose of this work is to investigate the dynamical behavior of a cooperative spin system, particularly, the 
dynamical critical behavior. It is well known that a great progress in the understanding of critical dynamics has been 
made, since Glauber [ l| and Kawasaki P| completed their pioneering work on the time-dependent Ising model. In the 
past two decades much research has been devoted to a better understanding of dynamical behavior of various systems, 
and many theoretical methods have been applied and developed |^|-|l9|. Among these studies, the main attention has 
been focused on the discrete spin systems, such as the Ising model and the Potts model, and a little on the continuous 
symmetry 0(n) spin systems. Nevertheless, as far as our knowledge goes, only a few analytical results were presented. 
In the present work and the previous paper pcfl , we are focusing on obtaining the exact analytical results. This is our 
I ' main motivation. 

The Gaussian model is a variation of the Ising model. It is a uniaxial continuous spin model that shows different 
■ static critical behavior from the Ising model. Although its static critical properties have been investigated clearly, 
little attention has been paid to dynamical critical behavior. This is also the reason we study the kinetic Gaussian 
model. Within the framework of Glauber dynamics in our previous paper pp| ], we have obtained dynamical critical 
exponent z = \ jv = 2 at the critical point K c — b/2d based on rigorous analytical derivation. 

To our knowledge, only the kinetic Ising model with time-dependent external field has been investigated in detail 
plj , p2| . The present work is attempting to investigate further the dynamic behavior of the kinetic Gaussian model 
with time-dependent external field. This paper is organized as follows: In Sec. |J, we fast summarize the basic theory 
of the single-spin transition critical dynamics, and then derive the fundamental equations of the kinetic Gaussian 



model in a periodic external field. In Sec. Ill, an isotropic d-dimensional hypercubic lattice Gaussian spin model 
I ■ is treated by Fourier's transformation. We exactly obtain the local magnetization and the equal-time spin-pair- 
correlation function. The critical characteristics of the dynamical relaxation r q , the complex susceptibility q), 
and the dynamical response of the system to the time-dependent external field are investigated. Finally, we end up 

O 



the paper with concluding remarks in Sec. IV 



^ I II. FUNDAMENTAL EQUATIONS 

U ■ 

A. Basic theory of the single-spin transition critical dynamics 

A single-spin transition critical dynamics based on Glauber's theory jjj, applying to both discrete-spin and 



continuous-spin systems, was presented in our previous paper |20|. For the sake of application here we only give 
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a summary. 

Spins system with Hamiltonian, TL ({cj}) , where <Xj is the spin of site i and can take discrete values or continuous 
values, interacts with a large heat bath with temperature T. The heat bath gives rise to spontaneous transition of 
spins via exchange of the energy. The probability of transition of the ith spin per unit time from one value <jj to 
another possible value &i is denoted by Wi{a% — > <7j). Under the assumption of single-spin transition, the probability 
distribution function P({<7i}, t) of the system, being in the configuration (ai, 02, • • • , cr/v) a t time t, is governed by the 
master equation 



(1) 



where the spin transition probability satisfies the following restrictive conditions: 
(1) Ergodicity: 



(2) Positivity: 

(3) Normalization: 



V(Xj,aj : Wjiaj -> a 6 ) ^ 0; 
Vctj,^ : Wj{<jj -> <7j) > 0; 



(2) 
(3) 
(4) 



(4) Detailed balance: 



in which 



Wj(aj ->■ 6-j) _ P eg (<Ti, ...,gAf) 

-> (Tj) P eg (<Ti,„. ,£7 3 -, ...,(Tjv) ; 



Pe q ({<?}) = ^exp[-pH(.{<r}% Z = £ exp[-/3H({a})], 



(5) 



(6) 



where P eq is the equilibrium Boltzmann distribution function, Z the partition function and 7i({er}) the Hamiltonian 
of the system. 

By use of both the master Eq. ([!]) and the normalized condition (Q), the time-evolving equations of the local 
magnetization and the equal-time spin-pair-correlation function can be expressed as, respectively 



where 



J t fa(*)) = - K 



dt 



(')) + E few* fa -**)) 

{o-j} \ (Tfc / 



((7fc(i)(^(i)) = -2 (a k (t)ai{t)) + 

+ <T i( t ) I E fa ~* * fe 



fa(*)> =£^(*)i'(fa} S *): 



^(fa},*), 



(7) 



(8) 



(9) 



(10) 
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Although the evolution starts with a certain initial state at t = 0, the system must be relaxed toward the final 
objective, which is an equilibrium state characterized by P eq = (1/Z) exp \—pH {{ex})] in the absence of time-dependent 
external field, via interaction with heat bath. In addition, it is usually considered that the transition probabilities 
of the individual spins depend merely on the momentary values of the neighboring spins as well as the influence of 
the heat bath. So, even if the transition probability cannot be derived exactly by means of microscope, the following 
form: 



W l (a l -> o-j) 



1 



exp 



-PH* 



(ii) 



is well chosen, where Qi is the coefficient determined by the normalized condition (|4]). Equation ([11]) means that 
the transition probability from <jj to &j only depends on the heat Boltzmann factor of the neighboring spins. If the 
system is in a periodic low-frequency external field, Eq. (fTy) is still a possible choice. 



B. The model and the fundamental equations 



To study further the dynamical behavior of the Gaussian spins system near the critical point, we put the system in 
a periodic low- frequency external field which may be regarded as electromagnetic wave p2| . The reduced Hamiltonian 
of the system under consideration is 



ffH = l\,,n,a , + 22 hi (t) (7i, 
(i,3) 1 



where 







K, 



Jij , ,,n Hi(t) H . 
-, hi \t ) = — = exp (iLut — iqxi) . 



kpT kpT 



(12) 



(13) 



the first sum goes over all nearest-neighbor pairs of lattice and the second over all sites. Unlike Ising spin model, the 
Gaussian model have two extensions p3| ]: first, the spin can take any real value in the range of (— oo, +oo); second, to 
prevent all spins from tending to infinity, the probability of finding a given spin between <7j and <Ji + d<7i is assumed 
to be the Gaussian-type distribution 



f(oi)d<Ji = \ — exp 
Ztt 



-2^ 



(14) 



where b is a distribution constant independent of temperature. 

In terms of those mentioned above, we can derive the fundamental equations of the kinetic Gaussian model in the 
external field. Following Eq. (|TT|), we choose the spin-transition probability as 



W l (a i -> Oj) 



1 



exp 



Ki.i+w&i&i+w + hi&i 



1 



exp [Ei&i] , 



where 



-w hi ) 



(15) 



(16) 



and means summation over nearest neighbors. As the spin variable takes continuous values, the summation for 
spin turns into the integration 



E 



f(a)da, 



then the normalized factor Qi can be determined as 

Qi = ^ exp (Ei&i) = J exp (Ei&i) f d&i = exp (l^J 



(17) 



(18) 



3 



and another useful combination formula can also be written as 



Substituting Eq. ( |T^ ) into the time-evolving Eqs. (0), (j|), one gets 

-(ai(t)) = -(ffi(t)) + ^^if M+ „,(cx i+ro (t)) + J/ii (t) , 



(19) 



(20) 



w 



+-[hi (t) {a (t))+h {t) (<n(t))]. 
Equations (EG) and (pll) are the fundamental equations of the kinetic Gaussian model in the external field. 



(21) 



III. EXACT SOLUTION 



For a d-dimensional isotropic hypercubic lattice, the dynamic equations of the local magnetization ( |20| ) and the 
spin-pair correlation function (pll) can be rewritten as 



(22) 



where 



and 



(^l(*)o-j(*)> = -2<oi(*)oi(t)) 



^^(n-nXoiCtJcrj^ + ^^Cri-rjXoiWoiCt)) 



[(a i (t))e- t -^- J + (a j (t)) e 



q =(91,92, ■ ■ • ,3d),r =(xi,x 2 , . . .,x d ),i = (ii,i 2 , ■ ■ .,«d),j = (j'i,j2, ■ • • , jd), 



j K, nearest-neighbor-pair, 
K( - ri ~ T i> ~ i 0, others. 



Introducing Fourier's transformation 



(M(q',i))=^(a i (i))e^' r S 



(23) 



(24) 



(25) 



(26) 



which satisfy 



(G(q', q", t)) = J] <a k (*) cr, (i)) e «^W. 

k,l 



P *(q-q )- r i — a , 

7VZ^ e -°q,q', 
l 



(27) 



(28) 



we have 



N ^ 
q 



l.k- 



1 - - b K{4) 



(29) 



(30) 



where 



eft 



1 - ±tf(q')) + (l - ^(q"! 



(G(q',q",i)> + 

[(M(q' ) t))5 q , q » + (M(q",f)) e^ 4 



(G(q',q",i)) 



K(q) = J2 K ( r i ~ rj)e ,;q (ri - r j ) 

i 

d d 

= KJ2 (e tq ' a + e-* q ' a ) = IK ^cos(gia), 



(31) 



(32) 



K(0) =2ifd, 



(33) 



and a is the lattice constant, qi is the ith component of the wave vector q, and d is the spatial dimensionality. 

First we solve Eq. (|3^). Obviously, it is a first-order linear inhomogeneous differential equation with the canonical 
form 



dy(t) 
dt 



P(t)y(t) =Q(t), 



and its general solution is 



where 



V (*) = ^77 



n(to)y(t )+ I KOQiZW 

to 



/i(t) = exp 



P(t)dt 



Applying it to Eq. (30), one can get the following exact solution 

t 



(M(q',i)) =exp 



<M(q', 0)) + J exp (J- + iwej 

[(M(q',0)) -ff x(^q , Mq,q']expf-— ) + #oX(w, qO-W^ 



where 



and 



1 - (1/6) K (q') i _ (2Jf/6) ^)« =1 cos(^a) 



X(w,q) = -7- 



6 l-(2X/6)EtiCOs(^a) 



(34) 



(35) 



(36) 



are the wave-number-dependent relaxation time and the frequency- and wave-number-dependent complex suscepti- 
bility, respectively. 
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From expression ( |3q ) we can see that r q < is finite for q' 7^ as the temperature approaches the critical point T c 
(K c = J/kpT c = b/2d), while it becomes to infinity for q' = 0. 

Substituting the solution ( |34| ) of the local magnetization into Eq. ( |3l| ) , the time evolution equation of the spin-pair 
correlation can be rewritten as 



(G(q',q",t)) + 



dt 

b 

NHoP 
b 



1 



1 - -K(q') + 1 - T X(q' 



1 



[<M(q', 0)) - H oX (u>, q')<5 q , q '] <J q , q » exp 



t 



(G(q',q",t)> 

+■ iuit 
t 



[<M (q", 0)> - H oX (u;, q")W] <W CX P — + iwt 
[fT x(w, q') + ffox(^, q")] -WW^'- 



(37) 



Equation ( |37j ) is also a first-order linear inhomogeneous differential equation. One can give its general solution 
(G(q',q",i)) 
= exp(-- tW( q ', q ", )) 



+ NH^l [(M(q // )0)) _ H x(o;,q")<5 q , q »] <W J* exp (Jj+™Z) # 



[fl0X(w. q') + ^0X(W, q")] <W<W 



b 
t 

x / exp 

= [(G(q',q",0)) - (M(q',0)) J ffox(w,q")W' " (W, 0)) H oX (fi>, q')W 

+^0X(w,q')x(^,q") ( 5q,q'^q"] eX P ( " ~ 



V r q' ' q' 

+ [(M(q', 0)> - H 0X (U, q')«q,q'] #0X^, q")*q,q« eX P 



+ [<Af (q", 0)> - i?ox(^, q")<Sq,q»] #ox(^, q')<W exp + iut 



r q " 



+Hlx{", q')x(w, q")<W<We 2M '- 
To make the solution an explicit one, we note that the factor exp (— t/V q ) can rewritten as 



exp 



exp 



1 — cos(fta) t 



e * TT exp 

i=i 



2K e l ^ a + e~ l ^ a 

T 1 2 



in which 



exp 



2Jf e^ Q + e" z « ia 
2 



is just the generating function for the first-kind imaginary argument Bessel function 

00 

exp^A + A- 1 )] = J2 A<,J «W. 



(38) 



(39) 



where I a (x) is the first kind imaginary argument Bessel function. Hence 

-p (-Pi =«-*n e (^r%(^ 



£=1 n^ — — oo 



i—l rti — — OG ^ ' 



for convenience, where the summations for nx, ■ ■ ■ ,nd from — oo to oo are denoted by 5Z n . Then Eqs. (34) and 
can be rewritten as 



(M(q',t)) = [(M(q',0)) - ff x(w, q')<W 
xe 



(G(q',q",t)) 

e' 2t [<G(q', q", 0)) - (M(q', 0)) ff oX ( W , q")<W' " < W> 0)) #ox(^, q')<W 



X 



ii . i ri 



+e- t e^ t [<M(q' ; 0)) - i? x(w, q')<W] H oX (u, q")<W 



x 



+e- t e- t [<M(q", 0)) - ff x(w, q")<Wl ffox(w, q')* q q ' 



b J d \ b 

II / \ 

+tf 2 X( W , q') X (w, q / ')*q,q'*q,q" ea<W *- 

Taking the inverse Fourier transformation 

<M*)) = ^E< M (^)> e ~ iq '' rk > 

q' 



K(*)<n (*)> = ^2 E <G(q',q",t))e-^'- r -^"- 
and using the following relation: 



iV 2 

q'.q" 



q 

the local magnetization and the pair correlation of the d-dimensional hypercubic system can be written as 



1 , s ,„ r (IK \ (IK 

■ -H 0X (uJ,q)e^I kl+ni ^—tj---I kd+nd i—t 

-Iff oX ( w ,q) e -*-^, 
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and 



(a k (t) ay (t)) 



]T (°) (°)> ^-n, ( ^ ) (~T' t l ■ ■ ■ J k«-n d ( ^« ) ^-m d ( ^ 



L- 2i J ffoX(^ ) q)e- <q -'^K (0)) 



Xhi- ni I — * I £ 



II . Ill 



t ■■■h d 



2K 



b 7 " c<i_nd V b 
-±e- 2t H oX (u, q)e- 4q - rk £ (a m (0)) 



2K 



T (IK \ /2if . 

1 V ~&~ / (l ~ m i I "ft -4 ' " 

+ ^ 2 %V(-,q)£e*' 



2K 



2K 



2K 



2K 



X/fe i+«i I — 1 J 7i i+"h ( — t ) ■■■ I k d +n d I -yt J Il d+md I — 1 



2^ 



2# 



2 

iV 2 "' 



c if x(^ q)e" iqri E ^ (°)> ^-"i (y*) ' ' ' ^-n d (^t) 



6 



-_tf V( w ,q) e ^-^ 



(47) 



respectively. Because n.j (or m,) can take any real value in the region (— oo, oo), the summations for rii and — n, (or 
for mi and —mi) are equal. In addition, the summation indexes n and m can exchange each other. Then, the Eqs. 
( flit ) an d (HI) can be rewritten as 



n 

i 



((7 n (0))-lff oX ( W) q)e^' 



2# 



2X 



iV 



#ox(^,q)e 



iujt— *q-r k 



(48) 



(<r n (0)a m (0)) + _tf o y( w , q )e-^ r «- 



2X \ /2if 



.'2/r 



2tf 

Id-ma \ 



E(^(0))e^/ fel _ ni 



2KT 



2K \ (2K 



2K 



-e-'e^Hoxiu, q)e-^ r . £ [(<r n (0)) - flbx(w, qK^] 



-f 



xI ki-n 1 \JY t ) '" Ik "- n d y 5 

+ ^iJ 2 X 2 ( W ,q)e 2 ^-^ r ^S 
where (cr n (0)) and (cr n (0) a m (0)) correspond to their initial values. 



(49) 
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So far we have obtained the exact solutions of the d-dimensional kinetic Gaussian model in the periodic external 
field. Letting Hq = 0, one can get 



<CT fcl ... fcd (t)) = 



OO 



/ / u ( 2K \ ( 2K \ 

/ J (c r „ 1 ...„ d (0j) Ikx-m [~ t J -"ha-m [~ f I ' 

m ...,rid.= — oo ^ * ^ 



(50) 



{(T kl ... kd {t)(j h ...i d {t)) = e 



xJ, 



E E 

ni j ...,n c j j — — oo mi ) ...,m t j=-oo 

,'2K\ _ /2if 

fci — n. 



(CTn 1 ...n d (0)(T mi ... mti (0)) 



ft -t 1 ifc-m, I -y*) ■■■h i -n d I /id_md IT* 



(51) 



Equations (£>0|) and ([51]) are just the exact solutions of the d-dimensional kinetic Gaussian model in zero external 
field, which agree with the results obtained in our previous paper IPO j. 

To understand the physical meaning of the Eqs. (^) and (|51|), it is interesting to investigate the long-time 
asymptotic behavior of the local magnetization and the spin-pair-correlation function. By use of the asymptotic 
expansion expression of the first kind of imaginary argument Bessel function 



I v (x) 



E 



(-)"(«, n) e -^+(«+2)" ~ ( v , n) 



2nx (2s) r 

n— v ' 



'2ttx 



E 



(2x)> 



, (— 7r/2 < arga; < 3n/2), \x\ — > oo, 



(52) 



where 



(v,ri) 



nW (| + v - n) ' 



one can get 



(*)) - E 



^n(o))--i?ox(^,q)e iq - r " 



-t/T 



N 



H x(u, q)e 



iut— *q-r k 



(53) 



where 



(tTk (t) 01 (t)) - XI 



(0) a m (0)) + jpHtfiw, q)e-^ r »- 



iV 



>q-ri- 



(^n(0)) 



-2t/r 



^0X( W ,q) e ^-' iq - r 'X [(^n (0)) -ff x(^,q) e - 4q - r "] 



(54) 



1 - 2fcd/6' 



(55) 



is the relaxation time of the system. From Eqs. (53) and (|54j ) we can see that both the local magnetization and the 
spin-pair-correlation function consist of two parts: one decays with t, and the other vibrates with t. Since r increases 
to infinity as the temperature approaches the static critical point T c (K c — J/kpT c = b/2d) , the decay term will occur 
critical slowing down phenomenon. 

We now turn on the response of the system to the time-dependent external field. According to the general theory of 
linear response, the complex susceptibility x(w, q) is expressed in terms of the equilibrium correlation of magnetization, 
namely |23] 



x(w,q) = x(o,q) 



P 1 Jo 



(M(-q,0)M(q,t)) f 



' iuJt dt, 



(56) 







where 



X(0,q) 



kpT 



(M(-q,0)M(q,0)\ 



and (•••)_ denotes the average over equilibrium distribution. 
Because 



N/3 
~b~ 



1 - 



IK 



Eqs. (|5q ) and ( |57[) , therefore, mean that 

(M(-q,0)M(q,t)) e = (M(-q, 0)M(q, 0)) e exp 



and 



(M(-q,0)M(q,0)) e = 



iV 



N 



6 [1 - ^ (q)] 5 



where 



l-i^(q) l-3fEtiCOstea) 



(57) 



(58) 



(59) 



(60) 



(61) 



It is interesting to note that as the temperature approaches the static critical point T c , for q = the static spatial 
correlation diverges, while for q ^ it remains finite. However, whether the singularity occurs or not, the dynamic 
responses strongly depend on the frequency w and the wave vector q of the external field. 



IV. CONCLUDING REMARKS 



In this paper, the single-spin transition Glauber dynamics is used to investigate the kinetic Gaussian model in a 
periodic external field. We have exactly obtained the local magnetization and equal-time pair correlation function of 
the cZ-dimensional isotropic hypercubic lattice Gaussian model by using Fourier's transformation. The related critical 
dynamics characteristics of the system are discussed. 

The master equation (Q) with the transition probability given by Eq. (|l"l| ) conserves the important features of a 
cooperative system. When the system is in a time-dependent external field, the dynamical model itself cannot be exact 
at very high-frequencies, it is only suitable for the case of low frequencies. In fact, the existence of a high frequency 
field weakens stochastic motion and makes a thermal equilibrium state with canonical distribution impossible. Even 
so, it would not bring any impact for the characteristic behavior of the system at low frequencies. The present work 
has clearly shown that the local magnetization, the equal-time pair correlation, and the dynamical responses of the 
system to a time-dependent external field strongly depend on the frequency and the wave vector of the external field, 
and when uj — > and q — * 0, they approach static results. 

The Gaussian model is certainly an idealization, but it is interesting and simple enough to obtain some fundamental 
knowledge of dynamical process in cooperative systems. Furthermore, although it is an extension of Ising model, the 
Gaussian model is quite different from Ising model in the properties of the phase transition . It is well known that 
in the equilibrium case the Gaussian model is exactly solvable on translational invariant lattices |pq| . Meantime, as 
we have done, the Gaussian model is also exactly solvable in dynamical case pcfl . Finally, we can anticipate that the 
kinetic Gaussian model will be a starting point to study the kinetic s 4 (or </> 4 ) model. 
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